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Abstract. 

Standard stochastic Loewner evolution (SLE) is driven by a continuous 
Brownian motion, which then produces a continuous fractal trace. If jumps arc 
added to the driving function, the trace branches. We consider a generalized SLE 
driven by a superposition of a Brownian motion and a stable Levy process. The 
situation is defined by the usual SLE parameter, k, as well as a which defines 
the shape of the stable Levy distribution. The resulting behavior is characterized 
by two descriptors: p, the probability that the trace self-intersects, and p, the 
probability that it will approach arbitrarily close to doing so. Using Dynkin's 
formula, these descriptors are shown to change qualitatively and singularly at 
critical values of n and a. It is reasonable to call such changes "phase transitions" . 
These transitions occur as k passes through four (a well-known result) and as a 
passes through one (a new result). Numerical simulations are then used to explore 
the associated touching and near-touching events. 



1. Introduction 

The scaling limit of many two-dimensional critical lattice models and growth models 
encountered in statistical physics may be described in terms of fractals. The 
boundaries of the Fortuin-Kastelyn clusters in critical g-state Potts models and critical 
percolation models, for example, in this limit are known to be fractal, conformally 
invariant curves EJE]. Stochastic Loewner evolution [212] (also described as Schramm- 
Loewner evolution and abbreviated as SLE) is a rigorous mathematical tool for 
producing and studying stochastic conformally invariant curves in the plane (see Rcfs. 
jSJ|l)l[7||H] for review). It was conjectured, therefore, that SLE is a description of such 
statistical systems. In several special cases this conjecture has been rigorously proven. 

SLE is based on the Loewner equation taken from complex analysis UJ. This 
equation contains a driving (or forcing) function, which determines all the 

properties of SLE. A few reasonable requirements constrain the choice of this function 
to a scaled Brownian motion, which leaves only one free parameter. One of these 
requirements (that the forcing function is continuous) ensures that the curves produced 
by SLE do not exhibit branching, thus leaving many systems, such as branching 
polymers, out of the picture. In this paper, we generalize SLE by dropping the demand 
that the forcing stochastic function be continuous, but keeping the requirement of 
stationary and statistically-independent increments. This leads to a much broader 
class of forcing processes including, in particular, the so called Levy processes 
^3^2 EH ESI- This generalization might be a useful description of tree-like stochastic 
growth. 
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Standard SLE has been studied in various geometries under the names of chordal, 
radial and other kinds of SLE's. In this paper we restrict ourselves to the chordal 
situation (growth in the upper half plane). We believe that our generalized SLE in a 
radial geometry may be relevant for description of the diffusion-limited aggregation, 
and plan to study it in a separate publication. 

In the next section, we define and describe the problem, in terms of the parameters 
which define the forcing. The forcing includes both a scaled Brownian motion and 
a stable Levy process. The parameters arc n, which sets the normalization of the 
Brownian term, c which sets the normalization of the Levy term, and a which 
determines its shape. This qualitative section ends with numerical simulations of 
the traces, the geometrical structures generated by the SLE. 

The final section of the paper gives the analysis of the short-distance properties of 
the traces showing both analytically and numerically that the traces have qualitative 
change in behavior as k and a each pass though critical values, respectively at four 
and one. The transition at k — 4 is quite analogous to a well known transition in the 
standard SLE [3]. The latter phase transition, at a = 1, is entirely new. This paper 
focuses upon the short-distance manifestations of this phase transitions in the chordal 
geometry. The accompanying large-scale features (which may be different for different 
geometries) will be the subject of a subsequent paper. 

We point out that most of the statements about the behavior of objects related 
to SLE and its generalizations are probabilistic, and typically appended by qualifiers 
"almost surely" or "with probability one" in mathematical literature. We will not 
use these qualifiers explicitly assuming that they should be applied wherever it is 
necessary. 

Some technical details of our calculations arc given in the Appendix. 



2. Description of problem 

2.1. Loewner evolution 

The so called chordal Loewner evolution is described by a differential equation obeyed 
by a function g(z,t). At each value of the "time" t, this function is a conformal map 
from a subset of the upper half of the complex z-plane (which we refer to as the 
physical plane) to the entire half plane H (mathematical plane). Specifically we write 

2 

dtg(z, t) = —-— -Tir, g(z, 0) = z. (1) 

g[z,t) -£(t) 

Here £(t) is the forcing, mentioned above. Alternatively, the function h(z,t) = 
9{z,t) -£(*) obeys 

d t h(z, t) = -A- - dt&t), h(z, 0) = z, (2) 
h(z, t) 

assuming that £ vanishes at t = 0. 

For any point z £ H these equations are valid up to the swallowing time t(z), 
which is defined as the first time when g(z,t) = £(£) or h(z,t) = so that the right 
hand sides of Eqs. (|TJ and become singular. The hull of the Loewner evolution is 
defined as the set of all points which become singular (are swallowed) in this way up 
to time t. At any given time, the z's that newly enter the hull may represent a single 
point or a more complicated set, sometimes including a subset H with a finite area. 
See Figure CI 
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Figure 1. The trace and the hull for a touching event. Here £(f) = 2\/6(l — t) 
for t 6 (0, 1) and zero elsewhere II 51 . Left shows the situation just before touching 
(t — » 1 — ); right shows the situation just after (t > 1). The trace is the thick dark 
line. The hull consists of that line plus the grey area. That area is added to the 
hull at t = 1. Note that there is a continuum of points added to the hull at the 
time of touching, but only one of these is on the hull's boundary and is the point 
7(1) of the trace. 

Another geometrical object considered in the SLE literature is the growing tip 
7(i) of the hull. The point z = j(t) is picked from among the (perhaps many) points 
which newly enter the hull at time t. Among these, this point is unique for having 
the property of being on the boundary of the hull at the time t, see Figure More 
formally, this point is defined by @] 

7(i) = lim/i _1 (z,i), (3) 

z—*0 

where the limit is taken within H. The union of the points j(t) up to a time T is the 
trace of the evolution up to T. 

The choice of the process £(t) defines the properties of the hull and the trace. 
In particular, if has a discontinuity (jumps), branching or discontinuity occurs 
in the trace and the hull may become disconnected. For a sufficiently smooth (a 
precise condition involves Holder continuity, see Refs. ^]E3E|) the trace is a curve 
parametrized by t. Moreover, if the trace is a simple curve, that is a curve without 
double points and without touchings of the boundary, it then coincides with the hull. 

2.2. The forcing 

Standard SLE |2j is the Loewner evolution with scaled Brownian forcing 

m = vnB(t), (4) 

where B(t) is a standard Brownian motion. In physical terms this means that the 
"time derivative" is a white noise with intensity k: = K $(t — 

The increments of in distinct time intervals are stationary and statistically 
independent. These properties of ensure that the process g{z,t) is conformally 
invariant in the sense that the time evolution is consistent with composition of 
conformal maps. In particular, the ensemble of g(t + s) can be obtained by functional 
composition of g(t) with g(s), each drawn from the correct ensemble. 
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If we want to generalize Eq. but maintain the functional composition property, 
we may only discard the requirement of continuity of the forcing £(t) but keep its 
increments to be stationary and statistically independent. This leads to a broad class 
of stochastic processes, including the so called Levy processes ^3^2 ED- a ls° called 
Levy flights in physics literature ^31 • The simplest of these are the so called stable 
Levy processes L a (t) characterized by a single real parameter a in the range < a < 2. 
The process L a {t) is composed of a succession of jumps of all sizes. In a time interval 
dt, L a (t) jumps by an amount between x and x + dx with the probability proportional 
to 

dtdx-. — 777— — . (5) 

More precisely, the probability distribution function of L a (t) is given by the Fourier 
transform 

P(x,t)=J°° e -^ e -*|fe| a g. (6) 

The case of a = 2, formally corresponding to a Brownian motion, is, in fact, quite 
different from a < 2. We find it interesting to combine a Brownian motion with a 
stable Levy process as the input to our generalized version of the SLE. Thus we choose 
as our forcing the function 

Z(t)=VZB(t)+cV a L a (t). (7) 

Every time a discontinuity of the driving force occurs the trace develops a branching 
point. With Levy flights, therefore, there are infinitely many branching points on all 
size-scales. 

The stable Levy processes are self-similar ^0)- For any positive number b the 
two processes, L a (t) and & 1//a Lo,(i/&), are "equal in distribution", that is, statistically 
identical, which we denote as 

L a (bt) = tf' a L a (t). (8) 

Thus, a typical trajectory of a Levy flight looks like many clusters separated by long 
jumps. Each cluster, if zoomed into, again looks like many clusters separated by long 
jumps, and so on |13| . 



2.3. Scales for SLE with Levy flights 

We choose the driving force £ (t) in Loewner evolution to be of the form : 

2 

dtg{z, t) = —— =— - , , g(z, 0) = z. (9) 

g(z,t) - y K,B(t) - c L ' a L a (t) 

If c = this process is scale- invariant \7\ . The run of the Loewner evolution up to a 
time bt is statistically the same as the run up to a time t, rescaled by a factor of y/b. 
The properties of the growing hull do not change with time except for general size 
rescaling. In this Brownian driving force is exceptional. 

The addition of Levy flights changes the situation: time becomes significant. 
Consider a time rescaling t — ► bt, where 6 is a positive number. We construct the 
rescaled conformal map 

g b (z,t) = b~ig(b^z,bt). (10) 
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It satisfies the initial condition g^{z, 0) = z, and using the self-similarity property © 
we observe that gi> obeys an equation of exactly the same form as Eq. except that 
the value of c is changed to 

Cb = cb^' 2 . (11) 

The strength of the Levy process flows under time-rescaling, making SLE not scale- 
invariant. (Many calculations in the usual SLE with pure Brownian motion are based 
on scale-invariance El El 113 so that they cannot immediately be extended to the 
present case.) 

This situation can be phrased in the language of the renormalization group (RG). 
The critical behavior of a system described by an SLE with purely Brownian forcing 
is determined by a fixed point of an RG. Levy flights make the system non-scale 
invariant, and their addition to the forcing can be thought of as a perturbation that 
is relevant at the hxed point. We do not know at present what such perturbations 
might be in terms of microscopic models of critical systems described by Brownian 
SLE. One possible candidate is a long range interaction of some kind. 

The flow of c has a fixed point at zero, when we recover the result that SLE with 
pure Brownian motion is scale-invariant. For c > 0, the ratio of K and the effective c 
changes with time. At time t = 1 it is just ci = c. At time t the effective ratio of the 
two terms in the forcing is 



,!/« c i/a 



-off 



K l/2 K l/2 

Therefore, crossover behavior can be observed at a time 



(12) 



l/(2-a) 



e = . (is) 

At small times t < f the evolution is dominated by the Brownian motion, whereas 
at large times t 3> t* it is dominated by Levy flights. (If k were very large or very 
small one would have to consider the independent effect of the three terms in the 
denominator of Eq. but for k of order one this is the only crossover that need be 
considered.) There is a similar crossover in space at a spatial scale x* = y/t*, with 
the Brownian motion dominating the small scale events. 



2.4- Numerical representation of SLE 

We approximate the integration of the SLE by introducing a discretized version, 
where each realization of £(i) is replaced by a piecewise constant function with jumps 
appropriately distributed. Let Xi, i = 1,2,..., be a sequence of real stochastic 
variables. We define a realization of the discretized driving force to be constant in the 
interval ((j — l)r, jr] (where r defines the mesh on the time variable) with £(0) = 
and 

3 

Z(t)^t3=J2 Xi for U-l>r<t<jT. (14) 

i=l 

For this forcing we can exactly solve for the trace. Within each time interval (ij-i, tj) 
we can solve the SLE forward: 

9(zJt) = y/\g(z, (j - l)r) - &] 2 + 4r + (15) 
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or backward: 



g(z, (j - l)r) = ^/[ff(z,jr)-^-] 2 -4r + 



(16) 



Equations (|15fl . i|16[l define an infinitesimal conformal map and its inverse for each 
time interval j: 



The trace can then be calculated numerically as an iteration process of infinitesimal 
conformal mappings starting from the condition g(j(t),i) = as follows |18l I19j : 



The value of ^ is determined by the variables Xj which are randomly drawn from the 
appropriate distribution according to the properties of the driving force £(t). 

In Figures[2HHlwe show typical SLE traces for different types of driving force up to 
a time T. For Figures |2H3 we have set c = and used increments drawn from Gaussian 
random variables in order to show SLE traces for the standard case of Brownian noise. 
The noise realization for these two pictures is the same, multiplied by an appropriate 
k. As is well known @j,[Zj, the trace does not touch itself for n ^ 4, but does so for 
k > 4. The figures show these contrasting behaviors. 

For Figures 01 we approximate SLE with Levy flights using the method in |2()j . 
As expected, these figures show branching structures produced by the discontinuous 
jumps. The figures suggest that a forcing consisting solely of discontinuous jumps 
does not produce any self-intersections. 

Figures EHU show the combined effects of jumps and Brownian motion. Exactly 
the same noise realizations have been used as in the previous four figures. The different 
sources of change in X (Brownian and Levy) have simply been added to one another. 

Comparison between Figures EH3 and [2H3] exemplifies the crossover behavior 
discussed in the previous section. Take for example Figure The trace consists of 
branches which look like the ones produced by the Loewner equation with constant 
forcing. Whenever a jump occurs, a new branch starts growing off some point on the 
trace or the real axis. If now we add Brownian noise to the driving force (FigurcsEHHJ) 
the trace at the largest scales makes branches like in the case with only Levy flights. 
However, if we zoom in (see insets) the trace appears similar to the SLE trace with 
only Brownian forcing (Figures [2H3- 

3. Analysis of short-distance behavior 

3.1. Phases of SLE 

It is well known 01 that the standard Brownian SLE process is qualitatively different 
in its behavior in the three different intervals < k ^ 4, 4<k<8, and 8 ^ n < oo. 
In the first interval the SLE trace forms a fractal that neither touches itself nor the 
real line (see Figure |2J). In the second interval the trace touches itself and the real 
line (see Figure 01 ■ For large enough time it will touch down arbitrarily far from the 
origin. Eventually it will surround and swallow up any point on the real line no matter 
how far it is from the origin. In the third, all points on the real line are surrounded 
sequentially, but also the trace fills in all the surrounded areas. This kind of behavior 




(17) 



(18) 



7j = lijr) = h o / 2 o . . . o fjfa). 



(19) 
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Figure 5. SLE trace for Levy distributed forcing a = 1.3, c = 4, 300000 steps, 
T = 10-'\ T = 3 




Figure 7. SLE trace with both Levy flights and Brownian forcing a = 1.3, c = 4, 
ft = 2, 300000 steps, r = 10" 5 , T = 3, t* = 0.31 




Figure 9. SLE trace with both Levy flights and Brownian forcing a = 
ft = 6, 300000 steps, r = 10~ 5 , T = 3, t* = 0.53 



7, c = 4, 
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is not discussed here. The different behaviors for different values of k are described as 
different phases of SLE, basically because the behavior in each interval is qualitatively 
different from that in the others. The difference is topological in character. It can be 
seen in both the small- and the large-scale features of the trace. Here we shall discuss 
phases of the model containing both Brownian and Levy forcing, but concentrate 
mostly on small-scale features of touching and near-touching. 

Notice that the events when the trace touches the real line and the previously 
existing parts of itself, are equivalent in the statistical sense, which follows from the 
stationary increments of the forcing and the composition property of conformal maps. 
Since the driving force £(t) is stationary and SLE is defined by a first-order differential 
equation, we have the statistical identity 

h(h(z,s),t) = h(z,s + t). (20) 

Consider the trace grown up to a time s + 1. The part of the trace that developed 
earlier (up to time s) can be mapped to the real axis by h(z, s). Then the event of the 
growing hull swallowing (or coming arbitrarily close to swallowing) a point somewhere 
on its previously grown boundary is statistically equivalent to the event of the hull 
swallowing (or almost swallowing) a point on the real axis, and the probabilities of such 
events should be the same. Therefore, to determine the probability of swallowing for 
points anywhere on the boundary of the hull, it is enough to determine the probability 
of swallowing only for the points on the real axis. This was initially argued for purely 
Brownian SLE, and remains valid here. 



3.2. Using Dynkin's formula 

The main tool that we will use to determine phase transitions in SLE is Dynkin's 
formula [211 12"2] . It can be written for a generic stochastic process (not necessarily 
Levy), but we will restrict it to Levy processes and stochastic integrals of them, that 
is, processes X(t) defined by a Langevin-like equation 

d t X t = F(X(t)) + a t (Levy process), (21) 

with any drift function F. Dynkin's formula gives the average value of any function 
/ of the stochastic process X(t) at an exit time T: 

</(X(T))) = /(.to) + (jAf(X(t))dt). (22) 

Here xq is the initial value of the process: X(0) = xq, and A is the generator of the 
process, that is, the operator acting on the argument of the function / and representing 
time differentiation. By choosing f(x) to be a zero mode of the generator: 

A f{x) = 0, (23) 

we get rid of the integral term. Note that we may solve (j23|) with any boundary 
condition: any non-constant solution is sufficient. In the Brownian case this equation 
is differential and easy to solve. Thus, Dynkin's formula produces many results about 
Brownian motion, such as recurrence, transience, etc. |21j . 
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3. 3. Brownian forcing 

As an example of using Dynkin's formula we first derive the phase transition at k — 4 
for purely Brownian SLE. This calculation is usually done in literature with the aid 
of martingales 0]. 

Let us set c = in Eq. {JJ. As a function of time, h(z, t) is a continuous stochastic 
process with the initial condition h(z, 0) = z. The event in which z is swallowed by 
the hull at a time T corresponds to the event h(z, T) = on the mathematical plane. 
We choose z to be on the real axis in the physical plane, then h(z, t) is real. 

Let us consider a real x > for concreteness. On the mathematical plane, we fix 
points a and b on the real axis so that 

0<a<x<6<oo. (24) 

Let T be the exit time from the interval [a, b]. Being continuous, the process h(x,t) 
can exit [a, b] either through a (with probability p a ), or through b (with probability 
1 —p a )- Clearly, h(x,T) is either a, or b. Then Dynkin's formula becomes 

p«/(o) + (l-p„)/(6) = /(ar) > (25) 

and allows us to find p a , the probability for h(x, t) to hit a before hitting b. Now if we 
take the limits a — > 0, b — > oo, it becomes the probability for h(x, t) to hit in a finite 
time, which is the probability for the point x on the physical plane to be swallowed 
by the hull: 

p=]imJSm MzM. (26) 
In general, the order of limits matters here. If it is reversed, 

p = um l im /(*) - f[ b \ (27) 

is the probability for h(z, t) to come arbitrarily close to 0, that is, for the hull to come 
arbitrarily close to the boundary. 

The generator for the process h(z, t) driven by the Brownian motion J3J| is [2] 

Ab = \d x + |e (28) 
A zero mode of this operator is 

' B W-\ log \x\ for K = 4. {ZJ) 

Substituting this into Eqs. (|26II27|) we find that the answer is independent of z, thus 
the probability for the hull to touch the boundary 

f for k < 4, , , 

Pc = = { 1 for K >4. ( 3 °) 

For k = 4 the order of limits is important and we find that 

/ for k < 4, 

Pc =° = { 1 for k> 4. ( 31 ) 

Note that in order to determine the phases, only the behavior of a zero mode of A 
at zero and at infinity is necessary. Examples of traces of standard SLE are shown in 
Figures [5] and |3J Touching is visible for the k — 6 and apparently does not occur for 
K = 2. 
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3.4- Dynkin's formula with Levy flights 

The case of the generalized SLE is different from the Brownian case: now h(z, t) has 
discontinuities. Therefore, at the exit time T the value of h(z,T) does not have to be 
a or b: the process can jump and overshoot the boundary of the interval. To avoid 
this difficulty, we define T to be the exit time from the set 

S= [-b, a) U (a, b], (32) 

and also restrict ourselves to zero modes that are even functions of x. 

The process h(z,t) can exit S by either hitting a (no overshooting possible), or 
by taking a value beyond —b or 6, in which case the overshooting is possible. But we 
consider the limit b — > oo. The behavior of zero modes, considered below, is simple at 
infinity: they either diverge or are constant. With these changes complete, we return 
to Eqs. (H3|23l. 

The generator for the process h(z,t) defined by Eqs. J21 0) is (see (TJJ) 
A = Ab + Al, where Ab is the same as in Eq. <|28[) . and 

AM = - 2r( _ a ; OT ^ + V) - /(*)) (33) 

This operator contains the principal value integral, which ensures convergence at small 
y for a < 2. The meaning of this term is most transparent in the Fourier space, where 
it is equivalent to the multiplication by — c|fc| a . 

Once again, we need to find a zero mode f(x), obeying Af = 0. These zero modes 
are studied in detail in the Appendix. Here we summarize the results. 

First of all, we find by direct calculation that the function 

/*<•> -{£&:• it\ (a*, 

is a zero mode of Aj,. Comparing this with Eq. (|29|l . we conclude that Jb{x) = Jl{x) 
if 1 — 4/k = a — 1 , or 

a = 2--. (35) 

K 

Along this line on the phase diagram we have an explicit expression for a zero mode 
of the full operator A. 

Secondly, for arbitrary a and k we can transform the zero mode equation, Af = 0, 
into Fourier space: 

^ [{c\k\ a + |fc 2 )/(fc)] - 2kf(k) - 0, (36) 

where f(k) is the Fourier transform of f{x). A solution of this equation is 

2c 



f (k) = \k\- a [i + -\kr a ) (2 - a,K . (37) 



Apart from a multiplicative constant and a possible delta function at k = 0, this 
solution is unique. Its asymptotic behavior is 

/<*)-{ \f: 2 ; k k Z7: ^ 

Note that Brownian motion determines the behavior at large k (small distances in 
x-space), whereas Levy flights determine the behavior at small k (large distances in 
x-space). 
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Figure 10. The phase diagram of the generalized SLE J5J. The dashed line is 
the hyperbola given by Eq. 1351 in the text. Along this line an explicit expression 
for the zero mode of the generator A is known, see Eq. 1311 . 



Analysis of the asymptotic behavior of the inverse Fourier transform of Eq. (|37fl 
gives the following: 

!|a;| 1_ « for K 7^ 4 and x — » 0, 

log | a; | for k = 4 and x — ► 0, /„„n 

Ixl"" 1 for a ^ 1 and x -> oo, ( ' 
log |ar| for a = 1 and x — > oo. 

Strictly speaking, we are able to obtain these asymptotics only for k > 2. The trouble 
is that the function that behaves as |x| 1_4 / K with k < 2 for small x would lead to 
a divergent integral in Eq. (J2HJ, and so cannot be a zero mode of A. Nevertheless, 
on physical grounds we believe that there is nothing special happening at n = 2, and 
that our results for phases that we state next, apply to all values k > 0. 

The use of the asymtotic behavior l|39|l in the formulas l|2t)l I27fl gives the phases 
of SLE driven by the noise Q: 

• isolated trees k ^ 4, a < 1: p = 0, p = 0. An example of the trace is shown 
in Figure El The Levy flights have broad distribution. Long jumps occur often. 
The trace can jump far from the existing trees and then a new isolated tree starts 
growing. The chances of returning to one of the old trees are small, so the old trees 
stop growing. We hope to make this point quantitatively in a future publication. 

• a forest of trees k^4, 1 ^ a < 2: p = 0, p — 1. An example of the trace is 
shown in Figured The Levy flights rarely produce very big jumps, so that the 
trace starts a new tree close to the old one, or does not leave the old tree at all. 
The hull looks like trees growing close to each other with their branches almost 
touching. 
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Figure 11. We calculate which points on the real axis plane have been swallowed 
by calculating H(x) = min{h(x,t);t < T}. There is a direct correspondence 
between the areas swallowed by the trace on the physical plane and the process 
h(x, t) on the mathematical plane. 



• a forest of bushes n > 4, 1 ^ a < 2: p = 1, p = 1. An example of the trace 
is shown in Figure |H1 The Brownian motion, which controls the small scales, is 
stronger now. The trees become thicker and different branches touch each other. 

• isolated bushes k > 4, a < 1: p and p are equal and finite, that is, the 
probability for a fixed point x on the real axis to be swallowed by the hull at a 
finite time is a function p(x) of the position of this point . This function is defined 
by any of the equations (|26II27|1 (the order of limits is unimportant in this phase), 
where f(x) is the inverse Fourier transform of (|37[) . p(x) is a symmetric function 
decaying from p(0) = 1 to p(oo) = with the asymptotic behavior: 

p{x) ~ M" -1 forx^oo, 1 -p{x) ~ for x -> 0.(40) 

An example of the trace is shown in Figure Since p{x) varies between and 1, 
parts of the real axis will not be swallowed, leaving gaps among trees. 

As long as c > 0, the first three phases do not depend on c at all. This is an 
expected property, because unlike k and a, c flows under the renormalization ljll|l and 
p and p are constants. A phase diagram of the generalized SLE is shown in Fig. 1101 

3.5. Probability p(x) of swallowing on the real axis 

Figures 02 and [S] seem to show that the trace swallows whole areas and touches 
itself. In the numerical approximation where the forcing consists of discontinuous 
steps the trace doesn't exactly touch itself, however it comes closer to touching as 
t — ► and whole areas thus appear to be swallowed. We can then try to determine 
the probabilities for swallowing by setting a small-distance cut-off and looking for 
limiting behavior as the cut-off goes to zero. 
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Figure 12. SLE with Levy Flights and Brownian component(ft = 7.5, a = 
0.5, c = 10). Main - Large x behavior: Probability p(x) of a point x on the real 
axis not to be swallowed until time T = 10 (r = 10~ 5 , 1000000 steps, t* = 0.09, 
H c = 3 10 — 3 , statistics for 7000 independent noise realizations). Inset - Small x 
behavior: Probability 1 -p(x) (r = 10~ 6 , 1000000 steps, H c = 6 10~ 4 , statistics 
for 3000 independent noise realizations). In both cases the dotted line is the 
theoretical prediction (Eq. 1391 *1 which is followed accurately by the numerics for 
four decades of x. Deviations for small x < 0.003 are due to the finite time step 
r of the noise realizations; the data points approach the theoretical prediction as 
we approach the scaling limit (not shown) 



Iii order to measure p(x) we consider the real stochastic process h(x,t) on the 
mathematical plane defined by Eq. (J2J). A point x on the axis is considered to be 
swallowed at time T if |/i(a;, t)| < H c for some t < T. Our method for determining 
H c is shown in Figure 1111 For one particular realization of the noise we calculate 
H(x) = mm{\h(x,i)\,t ^ T} for x on the real axis. We use the same realization to 
draw the trace on the physical plane. From the trace we can infer which points of the 
axis have been surrounded by the curve and thus have been swallowed. In the figure 
the boundaries of a swallowed area on the real axis are marked by the arrows; the 
dashed line corresponds to the cut off H c , we see that all points at which H{x) < H c 
appear to be swallowed in the physical plane. We repeat for several noise realizations 
to determine the value of the cut off accurately. The existence of the cut-off means 
we will be unable to distinguish between p and p, so that the case a > 1, K < 4 (the 
forest of trees) will be similar to a > 1, K > 4 (the forest of bushes). For all cases we 

have taken T > t* = ( — J 

Figure IT^l shows the results of the probability distribution p{x) of a point on the 
real axis being swallowed by the time T — 10 for a case with a < 1, n > 4 (isolated 
bushes). The numerical results closely follow the theoretical prediction of Eq. (|40|l . 

In general, we reproduce the theoretical form of p(x) rather well in this phase, but 
less well in the other phases. This is not surprising, it is the effect of finite time T in 
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our calculation. In the phases with forests, points are surrounded by a process which 
includes multiple returns to the same region of x. In the simulations, not enough 
time is allowed for this multiple return to work itself out. Furthermore, our use of 
the cut-off means we are unable to distinguish between p and p, so that the method 
makes the forest of trees look similar to the forest of bushes. Thus, we shall have to 
put together a more refined method of analysis and do much longer runs to make this 
distinction. 

Discussion and conclusions 

The generalization of stochastic Loewner evolution, suggested in this paper, has 
many new features. The hull acquires branching and the growth becomes a tree- 
like stochastic process. The phases of such hull are very different from the usual SLE. 
Although the driving force of SLE is still self-similar, the growing hull is not, the force 
of the Levy flights increases with time, being a relevant parameter under size rescaling. 

The non-stationary nature of the growth considered in this paper is related to a 
technical difficulty in its study. Many properties of the standard SLE are obtained 
using the so-called random time change which, in the case of the Brownian forcing, is 
consistent with the form of the Loewner equation Q . In the case of Levy processes 
a similar random time change exists, but changes the nature of the Loewner equation 
© and, therefore, cannot be used to extract useful information. 

In terms of the geometry of the generalized SLE driven by Levy processes, the non- 
stationarity means that many such properties as, for example, the fractal dimension 
of the trace cannot be globally defined. Instead, they may possibly be defined locally 
at each stage of the evolution. The question of how one can characterize the geometry 
of the generalized SLE will be addressed in a future publication. 
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Appendix 

In this Appendix we study zero modes of the operator A = Ab + Al, where 



Direct approach 

A zero mode of Ab is given in Eq. (|29[) . To find zero modes of Al, let us consider the 
integral 




(A.l) 




(A.2) 
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This integral converges at large y if < a, at small y for any and any a < 2 (because 
of the principal value), and at y ~ —x for > —1. Thus, we assume — 1 < < a < 2. 
By rescaling y —> \x\y we reduce the integral (|A.2fl (for nonzero x) to 

I ,/3-a r f\y+l\ l3 -l 



1+a 



-dy. 



(A.3) 



The last integral here can be evaluated. First, we split the integration interval into 
three intervals and change variables in each appropriately: 

-1 r -t fixi, |„, , ji/3 _ ^ 



\v\ 



-dy 



(y i) f3 i 



,1+a 



dy + 



1 (i-vY 



-dy 



,1+a 



dy. 



The integrals here can be found in Rcf. Ej] . Some of them contain hypergeometric 
functions that have to be transformed further by applying the Kummcr's connection 
formulas (see Sec. 2.9. in Chapter 2 of Ref. jUJ). Combining all the results, we 
obtain 

T(0+l)T(a-0) , T(0 + l)T(-a) , T(a- 0)T(-a) 



a,0 



T(a + 1) 



T(0-a + l) 



+ -e- a ( 2 F l (-a, -0; 1 - o;e) + 2 F 1 (-a, -0; 1 - a; -e) - 2). 
a 

Finally, in the limit e — > 0, the last term here vanishes, and we get 

_ T(0 + l)T(a-0) T(0+l)T(-a) T(a-0)T(-a) 

aJi ~ r(a + i) + r(/?-a + i) + r(-/3) ' (A ' 4j 

Note that this expression vanishes identically if = a — 1, thus showing that 
fh{x) = |x| Q_1 is a zero mode of A^. 

Let us now consider the limit of — ► 0. Specifically, we note that 



J a (x) = UmyS 1 I a ^{x) 



In 









l+a 



dj/ = J Q |x| Q , 



J a — 



In \y + 1 

l + Q 



dy. 



The last integral can be obtained form Eq 
Expanding gamma functions we get 

0- 1 I a ,p = -(iP(l -a)- i/>(a)) - T(a)T(-a) + 0(0). 
a 

Finally, using properties of the gamma and digamma functions, we get 
J a 



7r na 
— cot — . 

a 2 



(A.5) 

4j) by carefully taking the limit — » 0. 

(A.6) 

(A.7) 



This expression vanishes at a = 1 showing that Jl{x) = ln|x| is a zero mode of Ajj 
for a = 1. 

Thus, we have shown that 

Mx) = { \x\ a ~\ a^l, 
is a zero mode of Al- 



(A- 
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Next we note that Jb{ x ) = Il( x ) ifl — 4/« = a — 1 or a = 2 — 4/k. Along this 
line on the phase diagram we have an explicit expression for a zero mode of the full 
operator A, which behaves for small and large values of x exactly as claimed in the 
text of the paper. However, this is restricted to the range of k > 2 (because a must 
be positive). This is, actually, more generally true. Indeed, the integral I a ^(x) with 
= 1 — 4/re converges at y near — x only for k > 2. Thus, our claims in the paper 
about the phases can be justified strictly speaking only for k > 2. 



Fourier transform 

After the Fourier transform the zero mode equation Af = becomes a differential 
equation 

[ 'c\k\ a + ^k 2 )f(k)]-2kf(k)=0, (A.9) 



dk 



Let us consider k > 0. Introducing the new function 



K . 



g(k)=[ck a + -k')f(k), (A.10) 



we rewrite the differential equation as 



dg _ 2k 2fc 1 " Q 
9 ck<* + -k 2 c +f fc2_ " 

d\n(c+ -k 2 - a ). (A.ll) 



(2 - a)K 

This is immediately integrated to give 

g (k)=(c+^k 2 -«) J ^. (A.12) 

The case k < is treated similarly, and we get in the end an even function (as 
expected): 



K 



f(k) = \k\- a (c+'-\kr a ) {2 - a)K . (A.13) 

Apart from a multiplicative constant and a possible delta function at k = 0, this 
solution is unique. 

Performing the inverse Fourier transform, we obtain a representation of our zero 
mode: 

r°° dk 1 f°° 

/(*)=/ 7re lkx f{k) = - dkcos(kx)f(k). (A.14) 

This integral diverges at the lower limit for a > 1, but this is a spurious divergence 
that is easy to amend. Namely, we can formally subtract a (possibly infinite) constant 
from f(x) replacing it by 

1 r°° 2 f°° /kr\ ~ 

f(x) = - dk(cos(kx)-l)f(k) = - dfcsin 2 ( — )/(fc).(A.15) 

TT Jo TT JO V 2 7 

This integral converges at the lower limit for all values a < 2, which is all we need. 

Situation with convergence at the upper limit is trickier. By Chartier's criterion 
(see Ref. [23); the integral (|A.14(l converges at the upper limit only if n > 2. This is 
the same condition that was obtained in the end of the previous subsection where the 
Fourier transform was not used at all. So, this seems to be a more problematic issue. 



-l 
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We also note that after the subtraction performed in Eq. I|A.15|) . the convergence at 
the upper limit worsens. Now the integral is convergent only for k > 4. But this 
can be ignored. We simply perform the integration for n > 4 and then analytically 
continue to smaller values of k. Convergence also improves when we take derivative 
in either of Eqs. (JA.14I IA.15|) with respect to the strength of the Levy term c. We 
can then in principle obtain f{x) integrating in c. 

Another comment is in order here. If we tried to obtain the zero mode Jb{x) for 
the Brownian SLE case using the Fourier transform, we would get 

f(x) = — e * kx \k\i- 2 = - dkcos(kx)k^- 2 (A.16) 



„ 2tt 

instead of Eq. (|A.14(1 . The last integral has the same behavior at infinity as the one 
in Eq. (|A.14|I . But the zero mode of Ag given in Eq. (|29|) exist for all values 

of K. This is another reason for us to believe that the value n = 2 is not special, and 
that our results for the phases of the generalized SLE apply for all values of k. 
Consider now the integral 

f(x,K,a)=J dkk- a (l + k 2 - a y 2 - c ' ) * sin 2 (^J, (A.17) 

where we ignored the overall constant coming from pulling a power c out of the integral, 
as well as rescaling of k and x. This integral can be evaluated explicitly through 
hypergeometric functions for such values of a that 

2-a = — , (A.18) 

where p,q £ N are mutually prime, see Ref. (2H|, Eq. 2.5.2.3. The Mathematica 
software can also do the integral with these values of a. The simplest of these 
corresponds to p = 1, q = 2, that is, a = 1. The corresponding integral is 

f(x,K,l)= f dk k- 1 (l + kji- 1 sin 2 (^) 



16(k + 4) \ k 2 k 4 / 

For all values of a given by Eq. (|A. 18|) the structure of the answers is similar. 
They are given in terms of hypergeometric functions p F q . The asymptotic analysis of 
these expressions can be performed with the help of Ref. [2U , and in all cases gives 
the result given in Eq. I|39|l . 
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